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Abstract 

We report a measurement of the ratios of the decay rates of the , and 
B^ mesons into exclusive final states containing a, J/ip meson. The final states 
were selected from 19.6 pb~^ oipp collisions recorded by the Collider Detector 
at Fermilab. These data are interpreted to determine the h quark fragmen- 
tation fractions /„, fd and fs- We also determine the branching fractions 
for the decay modes B+ J/i)K+, B+ J/il;K*{892)+ , B° J/i^K^, 
B^ J/^K* (892)0 _^ J/^^(1020). We discuss the implications of 

these measurements to B meson decay models. 
PACS Numbers: 13.25.Hw, 13.87.Fh, 14.40.Nd 
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Typeset using REVT^ 



I. INTRODUCTION 



The bound states of bottom quarks provide a laboratory in which we can investigate 
the behavior of the strong force (Quantum Chromodynamics or QCD) and the electroweak 
interaction |Jl|]. The lowest-lying bound states are the pseudoscalar mesons {B^ , and 
5^) formed by one bottom anti-quark bound to one of the three lightest quarks (m, d and 
s, respectively). The branching fractions {B) of these mesons into final states consisting 
of only hadrons (the fully hadronic decays) have been studied theoretically and have been 
shown to yield insights into the interactions that take place between a quark and anti-quark 
pair at short distance scales 

Experimental studies of bottom meson hadronic decays have been limited by their rela- 
tively small branching fractions (typically 10~^ to 10~'^) and the difficulty of detecting the 
final states. The most precise branching fraction measurements have been made at e^e~ 
colliders, where the B^ and B^ mesons are pair-produced at threshold [|],^. Bottom hadrons 
are produced copiously in high-energy proton-antiproton collisions and so it is possible 
to measure the branching fractions of bottom mesons into those fully hadronic final states 
that have distinctive final state topologies. We report a study of the branching fractions of 
bottom mesons into final states consisting oi a, J/ip meson and a light quark meson, using 
the Collider Detector at Fermilab (CDF). The data set consists of 19.6 pb~^ of 1.8 TeV 
pp collisions produced by the Fermilab Tevatron Collider. This work extends an earlier 
analysis of the same dataset |jlO| by incorporating an additional final state, J/'?/'0(lO2O), 
estimating the fragmentation fractions of B hadrons and providing improved measurements 
of the branching fractions of B^ and B^ mesons. Throughout this paper, references to a 
specific decay mode imply the charge conjugate mode as well. 

We have focused our study on the bottom meson decay modes that yield a. J/ip me- 
son that subsequently decays to a fi^fi^ final state. This results in a signature that we 
can readily identify using the CDF trigger system and provides the necessary rejection of 
other background processes. We have measured the observed cross sections times branching 
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fractions for the channels 



B+ J/^K+, (1) 

^ J/^K°, (2) 

B+ J/^/sT* (892)+, (3) 

5° ^ J/^K*(892)°, (4) 

S° ^ J/^0(1O2O). (5) 

The observed cross section for the decay mode B~^ J/ipK^ can be decomposed into the 
form 

cTo^s = ct{vV ^ h)UB{B+ ^ J/^ir+)e^^ (6) 

and similar forms can be written for the other decays. Here, aijpp b) is the bottom 
anti-quark production cross section and /„ is the probability that the fragmentation of a 
b anti-quark will result in a B^ meson. In a similar way, we define fa and fs to be the 
probabilities of a 6 anti-quark to hadronize and form a B^ and B^ meson, respectively. 
We will refer to these probabilities as fragmentation fractions and explicitly include in this 
fraction contributions from decays of heavier B hadrons into final states containing a B~^, 
B^ or Bg meson. The expression B{B^ J/ipK^) represents the branching fraction for 
this decay mode of the B~^ meson and e^^ is the efficiency of detecting the J/iIjK~^ final 
state. 

The fragmentation fractions into the different B meson states are not well known. If one 
assumes that these fractions are independent of the fiavor and energy of the quark initiating 
the hadronization process, then measurements of strange meson production in light quark 
fragmentation provide the most accurate estimates of fs Only one measurement of fs 
in the bottom meson system exists [|T2| and it has large uncertainties. No measurements 
have been made of /„ or f^ in the bottom meson system, despite the importance of these in 
many B hadron branching fraction, lifetime and mixing studies [|r^ . 
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The bottom quark production cross section cr(pp b) is not precisely known; the best 
measurements to date have uncertainties of order ±20% 0. The efficiency e^^ depends 
on an understanding of the b quark production properties and detector acceptance. We 
therefore present our measurements in the form of ratios of branching fractions in order to 
avoid introducing additional uncertainties due to the b quark production cross section and the 
final state detection efficiency. We then compare these ratios directly with phenomenological 
predictions of the relative B meson branching fractions. We also use our data to estimate 
the fragmentation fractions /„, fd and fs- 

We have organized this report as follows. In Section II, we describe the data selection 
and procedures that we followed to reconstruct the five decay modes. We present in Section 
III a study of the relative sizes of resonant and non-resonant Kir and KK contributions to 
the B meson final states. We describe the procedure used to determine the acceptance and 
efficiency corrections for each decay mode in Section IV. In Section V, we present the results 
of this study and conclude in Section VI. 

II. DATA COLLECTION AND SELECTION 
A. The CDF Detector 

CDF is a multi-purpose detector designed to study high-energy pp collisions produced 
by the Fermilab Tevatron Collider. It surrounds the interaction point with three charged 
particle tracking detectors immersed in a 1.4 T solenoidal magnetic field. The tracking 
system is contained within a hermetic calorimeter system that measures the energy flow of 
charged and neutral particles. Charged particle detectors outside the calorimeter are used 
to identify muon candidates. The detector has a coordinate system with the z axis along 
the proton beam direction. The polar angle 6 is defined relative to the z axis, r is the radius 
from this axis and is the azimuthal angle. Pseudorapidity is defined as ?7 = — lntan(^/2). 

The innermost tracking device is a silicon microstrip detector (SVX) located in the 
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region between 3.0 and 7.9 cm in radius from the beam axis. This is followed by a set of 
time projection chambers (VTX) that measure charged particle trajectories out to a radius 
of 22 cm. An 84-layer drift chamber (CTC) measures the particle trajectories in the region 
between 30 and 130 cm in radius from the beam. This tracking system has high efficiency for 
detecting charged particles with momentum transverse to the beam Pt > 0.35 GeV/c and 
\r]\ < 1.1, and the CTC and SVX together measure charged particle transverse momenta 
with a precision of ap^ ~ [(0.0066Pt)^ + (0.0009P|)2] ^2 (^ith Pt in units of GeV/c). 

The muon detection system consists of 4 layers of planar drift chambers separated from 
the interaction point by ~ 5 interaction lengths of material. Additional 4 layers of cham- 
bers are located outside the magnet return yoke (corresponding to 4 interaction lengths of 
material) in the central pseudorapidity region |?7| < 0.7 to reduce the probability of misiden- 
tifying penetrating hadrons as muon candidates. An additional set of chambers is located 
in the pseudorapidity interval 0.7 < \ri\ < 1.0 to extend the acceptance of the muon system. 
The muon system is capable of detecting muons with Pt ^ 1.4 GeV/c in a pseudorapidity 
interval \ri\ < 1.0. These and other elements of the CDF detector are described in more 
detail elsewhere W^. 



B. The J/tp Selection 

We selected the J/ final state using a three-level trigger system that identified collisions 
with two muon candidates. The first level trigger required that there be two track candidates 
observed in the muon system. The level one trigger track efficiency rises from ~ 40% at 
Pt = 1.5 GeV/c to ~ 93% for muons with Pt > 3 GeV/c. The second level trigger requires 
the detection of a charged track in the CTC using the Central Fast Track processor (CFT), 
which performs a partial reconstruction of all charged tracks above a transverse momentum 
of ~ 2.5 GeV/c. The CTC track is required to match within 15° in (p of the muon candidate. 
The CFT efficiency rises from 40% at a muon Pt = 2.6 GeV/c to ~ 94% for Pt > 3.1 GeV/c. 
The third level trigger requires that two reconstructed CTC tracks match with two tracks 
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in the muon chambers and that the dimuon invariant mass be between 2.8 and 3.4 GeV/c^. 
The efficiency of the level three trigger requirement is (97 ±2)% for J/ip candidates. There 
are 2.06 x 10^ dimuon candidate events that passed the level three trigger requirements. 

These events were further selected to identify a clean sample of J/ip candidates. We 
required that each muon candidate have a CTC track candidate with > 1.4 GeV/c. 
This track, when extrapolated to the muon chambers, was required to match within 3 
standard deviations of the extrapolation and measurement uncertainties with a muon track 
in the transverse plane (r-0) and along the beam axis direction. The two muon candidates 
were required to have opposite charges. We performed a least-squares fit of the two muon 
candidate tracks under the constraint that the two tracks come from a common point (a 
vertex constraint). We required the probability of this fit to be greater than 0.01. These 
requirements resulted in a signal of (7.89 ± 0.08) x 10^ J/ip decays on a background of non- 
resonant dimuon candidate events. The dimuon invariant mass distribution for this sample 
is shown in Fig. |I], along with an estimate of the background determined using same-charge 
muon candidate pairs. We performed an additional fit to the dimuon system, applying a 
vertex constraint and requiring that the dimuon invariant mass equal the world average J/ ijj 
mass of 3.09688 GeV/c^ [O. The confidence level of this vertex-plus-mass constrained fit 
was required to be greater than 0.01. 



C. Reconstruction of Exclusive Decays 

1. The B+ J/2jjK+ Channel 

We reconstructed the exclusive decay modes listed in Eqs. Ilj-^by forming charged particle 
combinations with the J/ip candidate. For the decay channel J/ipK^ , we considered 

every charged particle with Pj- > 1.5 GeV/c as a candidate and required the resulting 
candidate to have Pt > 8.0 GeV/c. A least-squares fit was performed on the three 
charged tracks forming the J/ ipK^ candidate by constraining the three tracks to come from 
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a common vertex, the invariant mass of the dimuon system to equal the world average J/iIj 
mass, and the flight path of the candidate to be parallel to its momentum vector in the 
transverse plane (a 2-dimensional pointing constraint). The confidence level of this least- 
squares fit had to exceed 0.01. We required the fitted transverse momenta of the muon 
candidate with lowest and highest Pt to be greater than 1.8 and 2.5 GeV/c, respectively. 
This ensured that the muon candidates were likely to pass the dimuon trigger requirements. 
In order to reduce the backgrounds from prompt J / ip production, we required the meson 
candidate fiight path to be pointing in the same hemisphere as its momentum vector (in 
effect requiring the B^ candidate's observed proper decay length, cr, to be positive). The 
interaction vertex position was determined by averaging the measured beam position over a 
large number of collisions recorded under identical Tevatron Collider operating conditions. 
The J/ipK^ invariant mass distribution is shown in Fig. 0(a). We have performed a binned 
maximum likelihood fit of this data to a Gaussian lineshape and a linear background term 
and estimate a B^ signal of 154 ± 19 events. The width of the signal was not constrained 
in the fit and resulted in a fitted mass resolution of 0.015 ± 0.002 GeV/c^, consistent with 
our expected detector resolution. 

The reconstruction of the other four decay modes was performed with similar criteria in 
order to reduce the systematic uncertainties resulting from the kinematics of the produced B 
mesons and selection biases. Identical requirements were made on the quality and transverse 
momenta of the muon candidates, constraints on the fits to the B decay topologies, B 
meson lifetimes and lifetimes. We allowed for small variations in the B and light quark 
meson transverse momentum requirements to optimize the expected significance for each 
channel. The significance is defined as Ng/ \/Ns + Nb, where A^^ is the expected number of 
events determined using a Monte Carlo calculation for a given integrated luminosity, and 
A*";, is the extrapolated background rate under the signal region using the observed B meson 
sideband background levels. This resulted in only modest differences in the Pt requirements 
from channel to channel, and did not introduce significant systematic uncertainties in our 
estimation of the B meson detection efficiency. 
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2. The 5° ^ Jj^K'^ Channel 

The decay mode J/ipK^ was reconstructed by searching for — * tt+tt" 

candidates using all pairs of oppositely charged particles. The daughter pions were required 
to have Pt > 0.4 GeV/c. To fully reconstruct the J/ipK^ decay, we performed a 

least-squares fit to the two pion candidate tracks and two muon candidates, constraining 
each track pair to come from common points, requiring the momentum vector of the 
candidate to point along its fiight path (a vertex and pointing constraint) and placing a 
J/ip mass constraint on the dimuon system. We also imposed a mass constraint on the 
dipion system, constraining the invariant tt+tt^ mass to the world average mass of 
0.4977 GeV/c^. The confidence level of the fit had to exceed 0.01. To improve the signal- 
to-noise ratio, we required the proper decay length of the candidate to be larger than 
0.1 cm and its transverse momentum to be greater than 1.5 GeV/c. The tt'^tt~ invariant mass 
distribution for the candidates that satisfy these requirements is illustrated in Fig. H(a), 
and shows a signal of (2.56 ± 0.05) x 10^ decays above a large combinatorial background 
(for illustration, no mass constraints were imposed in the least-squares fit to the charged 
tracks in this figure). To identify a clean B^ candidate sample, we required the J/ipK^ 
candidates to have Pt > 6.0 GeV/c and the candidates to have a cr greater than zero to 
reduce the combinatorial backgrounds from prompt J/ip production. The invariant J/ipK^ 
mass distribution for these candidates is shown in Fig. ^(b). A fit of this distribution to a 
Gaussian lineshape and linear background results in a total signal of 36.9 ± 7.3 B^ decays. 

3. The B+ ^ J/ ijK* (892)+ Channel 

We searched for the decay mode B^ J /ipK* (892)^ J/ipK^n^ by selecting a sample 
of J/ip candidate events containing a candidate. The criteria used to identify can- 
didates for the J/ipK^ final state were also used for this decay mode. We required the 
candidates to have Pt > 2.0 GeV/c. We considered all other charged tracks with Pt > 0.4 
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GeV/ c as vr"*" candidates and we combined these with the candidates to form all possible 
K*{892)^ K^TT~^ candidates. The combinatorial backgrounds to the i^'*(892)"'' decay are 
large, as illustrated in the KgTT~^ invariant mass distribution presented in Fig. ^(c). In order 
to identify a 5+ candidate sample, a least-squares fit similar to that imposed on the J/ ipK^ 
candidates was performed. We required that the confidence level of this fit be greater than 
0.01, and that the cr of the candidate be greater than 0.1 cm. The transverse momentum 
of the candidate had to exceed 6.0 GeV/ c and its cr had to be positive. In order to isolate 
a i^*(892)^ resonance, we required that the invariant Kgir^ mass be within 0.08 GeV/c^ of 
the world average K*{892)+ mass (0.8916 GeV/c^) This resuhs in the J/tpK^ir-^ invari- 
ant mass distribution shown in Fig. 0(c). A fit of this distribution to a Gaussian lineshape 
and linear background results in an estimated signal of 12.9 ± 4.3 decays. 



4. The B° J/tljK*{892)^ Channel 

Our data selection to reconstruct the decay J /ipK* [892)'^ —>■ J/ipK^n' proceeded 

in a similar manner. We formed combinations of all oppositely-charged track pairs, and fit 
the four charged tracks requiring that they come from a common decay point, constraining 
the invariant dimuon mass to the world average J/ip mass, and requiring that the flight 
path of the candidate be parallel to its momentum vector in the transverse plane. The 
confidence level of this fit had to be greater than 0.01 and the B^ candidate cr had to 
be positive. The combinatorial backgrounds to the K*{892)'^ decay are also large. This is 
illustrated in Fig. 0(d), where we show the K^tt^ invariant mass distribution for events that 
have transverse momentum of the K^tt^ system greater than 2.0 GeV/c. We defined the i?° 
candidate sample by requiring the Pt of the K'^n' system to be greater than 2.0 GeV/c and 
the resulting J/il)K^Tt~ system to have Pt > 8.0 GeV/c. We required the K~^7c~ invariant 
mass to be within 0.08 GeV/c^ of the world average K*{892f mass (0.8961 GeV/c^). The 
resulting J /ipK^TT~ invariant mass distribution is shown in Fig. |^(d). 

This peak also has contributions from K*{892)^ decays where the incorrect kaon and 
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pion mass assignments yield an invariant K^ti~ mass within the i^*(892)'' mass window of 
±0.08 GeV/c^. We used a Monte Carlo calculation, described in Section IV, to determine 
the relative fraction of such combinations and the shape of the resulting J/tpK^n^ invariant 
mass distribution. The signal shape was parametrized by two Gaussian distributions with 
the relative width, normalization and position of the second distribution determined by a fit 
to the J/ipK^TT~ invariant mass distribution predicted by the Monte Carlo calculation. The 
width of the second Gaussian was fixed to 3.3 times the width of the first, the normalization 
of the second was fixed to 0.08 times that of the first and the mean of the second Gaussian 
distribution was offset lower in mass by 0.0023 GeV/c^ relative to the mean of the first. This 
shape was then used in a fit to the observed J /ipK* (892)^ invariant mass distribution to 
determine the number of decays in our data. This procedure yields a signal of 95.5 ±14.3 
B° J/^ir*(892)° decays. 

5. The B° J/'0(/)(lO2O) Channel 

The search for the decay mode —>■ J/?/'0(lO2O) J/ipK^K^ was performed by 
considering as 0(1020) — * K^K'^ candidates all oppositely-charged track pairs. A least- 
squares fit of the ii~^^~K^K~ candidate system was performed, constraining all four tracks 
to come from the same vertex, constraining the dimuon invariant mass to the world average 
J/ ip mass, and imposing a 2-dimensional pointing constraint on the B^ decay. The confidence 
level of this fit had to exceed 0.01. A 0(1020) signal of (4.1±0.4) x 10'^ events is evident in this 
sample, as illustrated in Fig. ^(b). The combinatorial background was reduced by requiring 
the K+K- system to have Pt > 2.0 GeV/c, the J/i^K+K' system to have Pt > 6.0 GeV/c 
and the cr of the B^ candidate system to be positive. We defined our 0(1020) K^K~ 
candidate sample by requiring the K^K~ invariant mass to be within 0.0100 GeV/c^ of the 
world average mass (1.0194 GeV/c^) |T3|. This resulted in a sample with the J/iIjK~^K~ 
invariant mass distribution shown in Fig. 0(e). A 5° signal is evident on a relatively small 
background. A fit of this distribution to a Gaussian lineshape and linear background results 
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in a total signal of 29.4 ± 6.2 events. 



III. RESONANT AND NON-RESONANT DECAYS 

Clear signals for B meson production and decay are observed in all five channels. In the 
case of the three channels involving a K*{892) or 0(1020) resonance in the final state, the 
estimated number of B candidate events includes resonant and non-resonant contributions 
in the final state. We searched for evidence of a non-resonant Kn or KK contribution 
to the B meson signals by placing invariant mass cuts on the B candidate, removing the 
invariant mass cuts on the two-meson systems and examining the Kg7i~ , K~^tt~ and K^K~ 
invariant mass distributions. In order to account for non-i? background in the two-body mass 
distributions, we defined B mass sideband regions for the three samples, normalized to the 
estimated number of non-i? events as determined from the B invariant mass distributions. 
The signal and sideband regions are described in Table |I[ By allowing for a non-resonant 
contribution to the B decay rate, we will directly estimate the rate of resonant decays 
without having to assume that the rate of non-resonant decays is negligible. 

The K^TT'^ , K^TT^ and K^K^ invariant mass distributions are illustrated in Fig. ^ for 
the B signal and sideband regions (the shaded distributions are from candidates in the B 
sideband regions normalized to the background under the B meson peak in the signal re- 
gion). One sees from these distributions resonant signals for the fr*(892) and 0(1020) with 
no significant non-resonant contributions above the non-i? backgrounds. We quantified the 
amount of resonance production associated with the B signals by performing binned maxi- 
mum likelihood fits of the two-body invariant mass distributions to Breit-Wigner lineshapes 
convoluted with detector resolution, using the normalized B sideband distributions to model 
the shape and size of the background under the two-body resonance signals. The resulting 
numbers of observed events, Nsb-, are listed in Table |I[ As a cross-check, we also estimated the 
number of signal events by fitting the resonance signals to Breit-Wigner lineshapes convo- 
luted with detector resolution and background shapes described by second-order polynomial 
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functions. The resulting event rates, Nju, are listed in Table | and are consistent with A^^fe- 
Under the assumption that there are no non-resonant decays, we can also estimate the 
strength of the two-body resonant decay by correcting the observed B rates determined 
from the fits to the J/iIjKti and J/ipKK invariant mass distributions in Fig. 0(c)-(e) for 
the loss in efficiency due to the Kir and KK mass cuts. The presence of non- resonant 
Ktx or KK decays would result in corrected B decay rates systematically larger than those 
determined by Ngb or Nju. The mass cut efficiencies have been estimated using a Monte 
Carlo calculation to be 0.80 and 0.86 for the Kit and KK mass window cuts, respectively. 
The resulting B meson rates, Nwin-, are listed in Table |. We see no significant difference in 
the rates estimated by these three methods. We therefore conclude that we do not observe 
a significant non-resonant B — > J/ipKir or B^ J/ipKK decay mode. 

For the subsequent analysis, we choose Nsb as the best estimate of the rate of resonant 
production as it is least biased by potential contributions from non-resonant production. In 
addition, we have investigated the possibility that kinematic refiections of other B hadron 
decay modes could enhance our observed event yields, and have excluded such contributions. 



IV. EFFICIENCY CORRECTIONS 

We estimate the relative reconstruction efficiency for each B meson decay mode to convert 

the observed number of B events into ratios of branching fractions. We write the efficiencies 
for reconstructing B mesons as 

e - ej/V X ^geom ^ ^cr ^ {() 

= ej/v X efeL x ^cr x (8) 

e^*^ = ej/v, X ef,;^ x e^*^ x x (9) 



e = e 



J/V X e^eom X e,, X ei^ X (10) 
=J/V X ej,„„ X et X e<^, (11) 
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where we show the common contributions. The quantity ej/^ is the efficiency for triggering 
and reconstructing the J/ip ^ A*^/^^ decay. It is common to all decay modes. This also 
includes the combined efficiencies of the vertex and vertex-plus-mass constrained ffis of 
0.952 ± 0.006, which cancels out in our subsequent analysis. The quantities egeom are the 
geometrical efficiencies for finding the daughter mesons in the tracking fiducial volume, 
having the decay exceed the minimum Pt requirements on the meson and B systems given 
& J/ip candidate, and having the B candidate satisfy the constrained ffi requirements. The 
quantities ecr are the efficiencies of the proper decay length requirement on the B candidate 
in the different decay modes. The quantities e^, ex^, Ctt and are the efficiencies for 
reconstructing the , K^, vr" and mesons using the charged track information. 

In addition to these efficiencies, we correct the observed event rates for the relevant 
branching fractions into intermediate states, B{K^ tt+tt^) = 0.3430 ± 0.0014 and 

S(0(1O2O) K+K-) = 0.491 ± 0.009, both taken from Ref. [|13[, and the isospin weighting 
factors B{K* {892)^ K+n-) = B{K* (892)+ K%+) = 2/3. 

Since we are only interested in ratios of efficiencies, the common terms in these efficiencies 
cancel, reducing the overall uncertainties. These include the term ej/^, and the reconstruc- 
tion efficiencies ex or e^^ when they appear in both the numerator and denominator of the 
ratio. A number of other quantities do not necessarily cancel when calculating the ratio of 
branching fractions. In order to evaluate the relative efficiencies, we employed a B meson 
Monte Carlo calculation. B mesons were generated with a Pt spectrum predicted by a 



next-to-leading order QCD calculation [|15| using the MRS DO parton distribution functions 
|T6| . The b quark Pt was required to be > 5 GeV/c, and the b quark fragmentation into a 
B meson was modeled using the Peterson parametrization with the parameter e chosen to 



be 0.006 |T^. The B mesons were decayed using the CLEO B decay model [0 and a full 
simulation was used to model the response of the CDF detector, including effects due to 
the underlying event. The resulting Monte Carlo events were then processed with the same 
algorithms used to reconstruct the data. We used the reconstructed Monte Carlo events to 
estimate the geometrical acceptances egeom- These efficiencies are listed in Table ||. 
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The geometrical efficiencies include the effect of the B meson vertex and mass constrained 
fits. We determined that the efficiency of the ffiting procedure and subsequent confidence 
level requirements were independent of B decay mode by measuring the relative loss in sig- 
nal events when different event topologies were fit employing both a vertex constraint and 
a vertex-plus-mass constraint. We assigned a systematic uncertainty of 1% in the relative 
acceptances to this effect, which we estimated by comparing the relative loss of signal events 
in the different decay topologies. We also investigated the uncertainties associated with 
the model of the detector used to measure the reconstruction efficiencies. We verified that 
the detector simulations accurately described the interaction vertex distributions and de- 
tector geometry. We then compared the efficiencies determined using the complete detector 
simulation to those determined using a parametrized model of the detector. Based on this 
comparison, we assigned a 5% systematic uncertainty on the relative geometrical efficiencies 
to take into account any remaining uncertainties in the detector model. 

In principle, the B meson cr requirements have different efficiencies for each channel 
because of different momentum and vertex resolutions for the final states and possible dif- 
ferences in the lifetimes of the B meson states. We evaluated the efficiencies using the 
world average values for the lifetimes of the three B mesons [|I^] , using tracking detector res- 
olutions observed in the events in the B sideband regions. The results are listed in Table f|. 
We repeated the efficiency calculation varying the lifetimes by one standard deviation. The 
resulting variations in the relative efficiencies were assigned as systematic uncertainties. 

The meson reconstruction efficiencies take into account the charged track reconstruc- 
tion efficiencies, the efficiencies of the K'^ lifetime cut and the additional constrained ffis 
performed when a candidate is in the final state. The track reconstruction efficiencies 
were determined using both a full detector simulation and by embedding simulated tracks 
in real interactions containing J/ip candidates. The systematic uncertainties in the track 
reconstruction efficiencies were determined by varying the embedding techniques. The loss 
of K"^ mesons due to decays-in-flight was estimated using the full detector simulation. Be- 
tween 4% and 6% of mesons (depending on the B meson decay mode) decay within the 
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volume of the CTC, of which approximately 40% are correctly reconstructed. We assigned 
a 3% systematic uncertainty on ex due to the decays-in-fiight correction. 

The lifetime cut efficiency and the efficiencies of the vertex and mass constrained 
fits were determined by measuring the loss of signal events in both the proper lifetime 
distribution and the ti^ti^ invariant mass distributions. The efficiency of the lifetime cut 
was determined to be 0.958 ±0.007. The uncertainty represents the difference in efficiencies 
determined by estimating the loss of real decays using the tt^tt^ invariant mass and the 
cr distributions for the candidate samples. The fractions of candidates that satisfied 
the confidence level requirements on the vertex and vertex-plus-mass constrained fits were 
0.938 ± 0.017 and 0.983 ± 0.006, respectively. 

A number of additional checks were made to verify that correlations in efficiencies were 
properly taken into account. The variation of the dimuon trigger acceptance for the different 
B meson final states was determined using a Monte Carlo calculation that simulated both 
the detector response and the effect of the trigger. This resulted in a negligible uncertainty 
in the ratio of acceptances. Variations in the Pt spectrum of the produced B mesons 
could also result in a change in the relative acceptance. This was measured by varying the 
renormalization scale and the h quark mass in the Monte Carlo calculation of this spectrum. 
We assigned a systematic uncertainty on the relative acceptance due to this effect that varies 
from 1 to 5%, depending on the pair of final states being compared. The polarization of the 
vector mesons in the final state also has an effect on the relative acceptance. We varied the 
longitudinal polarization of the K*(892) meson in the B meson rest frame by ±0.10 around 
a nominal value of 0.75 and the 0(1020) meson longitudinal polarization by ±0.25 around 



the nominal value of 0.50 [|T9|. We assigned the resulting 2.5% change in acceptance as the 
systematic uncertainty due to this effect. 

We expect 5^, B^ and mesons to be produced both directly and through the pro- 
duction of excited B meson states that decay to the pseudoscalar mesons we observe. We 
investigated the effect such resonant production would have on the relative ratio of efficien- 
cies of the decay modes studied by performing a Monte Carlo calculation using the PYTHIA 
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program which models the production and decay of higher mass B meson states. In 



this calculation, we assumed that the relative production of B mesons with orbital angular 
momentum L and spin S was in the ratio 0.30 : 0.53 : 0.17 for L = 1 and 5' = 0orl: L = 
and 5 = 1: L = and 5* = [^. The change in the ratio of acceptances, relative to the case 



where only pseudoscalar meson production was assumed, varied from 1 to 4%, depending 
on the decay mode considered. We included this as an additional systematic uncertainty on 
the acceptance. 

The systematic uncertainties assigned to the relative efficiencies are summarized in Ta- 
ble |ITT[ These were combined in quadrature to determine the total systematic uncertainty 
on the relative acceptance for each pair of decay modes used in this study. 

V. RESULTS 

We present our results as a matrix of ratios of acceptance-corrected rates of B meson 
decays into the five channels. The observed numbers of signal events, listed in Table f|, 
were corrected by the detection efficiency for each decay. When we form the ten possible 
ratios of these acceptance-corrected event rates, the b quark production cross section and 
the common efficiencies cancel. The results are listed in Table Three of these ratios 



have also been determined using a different technique based on the same data set ||T0|- The 
values determined here are in good agreement with these previous results. Note that the 
two measurements of these three ratios are not statistically independent. 

The measured quantities are the ratios of the product of b quark fragmentation fractions 
and the B meson branching fractions into the specific final state. Thus our measurements 
can be written as 

- Tu Bim^jJW^) - ^-^^ ^ ^ ^^^^ 
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where the first and second uncertainties are the statistical and systematic uncertainties, 
respectively (henceforth the first and second uncertainties in measured values will represent 
the statistical and systematic uncertainties, respectively). 

We can use these data to constrain both the fragmentation fractions and the meson 
branching fractions. To extract the branching fractions, we will have to assume certain ratios 
of fragmentation fractions. Correspondingly, we will use phenomenological and theoretical 
predictions for the ratios of branching fractions to extract the fragmentation fractions. 

A. Branching Fractions 

1. The Bg Branching Fraction 

The ratios of branching fractions that involve the meson can be used with the world 
average values for the B'^ and B^ meson branching fractions into the four other final states 
13| to estimate the product of the ratio of fragmentation fractions, fs/{fu,fd), times the 



branching fraction B{B'^ J/ipcf)). The world average B^ and B^ branching fractions have 
been measured assuming that the fragmentation fractions and fd are equal for h quarks 
produced in T(4S') decays, and so this assumption is implicit in this calculation. 
All ratios that involve the B^ can be rewritten in the form given by the example 
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^ = fuB{B^ ^ J/^K+)Ri^, (22) 

which gives us four different measures of the ratio of fragmentation fractions and the 
branching fraction, using the world average values for the branching fractions |13| on the 



right-hand side of Eq. With the assumption of equal and B^ fragmentation fractions, 
we form the weighted average of the four estimates to obtain 

B{B^^ J/i)(t)) = (0.37 ± 0.11 ± 0.04) x 10^1 (23) 

\Jui Id) 

In order to extract B{B'^ — > J/ipcj)), we assume /„ = fd and use the value = 
(0.40 ± 0.06) fu- This value of fs represents the central value of the range of reported 
fs measurements [^,|T^, and the uncertainty has been chosen to cover half of the difference 
between the minimum and maximum values. It is also consistent with the suppression of 



strange hadrons observed in the production of light quark hadrons |2^. With these values 
for the fragmentation fractions, we determine 

i3(5° ^ J/^0) = (0.93 ± 0.28 ± 0.10 ± 0.14) x 10^1 (24) 

The first uncertainty is statistical, the second accounts for the systematic uncertainties 
associated with the ratio of branching fraction measurements and the third is the uncertainty 
associated with the value we have taken for fs- 

2. The and Branching Fractions 

We can also use these data to estimate the branching fractions E{B'^ —>■ J/tpK^), 
B{B^ Jli)K^), B{B+ J/i)K*+) and B{B^ J/^jK*^) using the world average val- 
ues for the branching fractions, our ratios of branching fractions and the assumption that 
fu = fd- For example, for the decay B^ — > J/tpK*^, we have three separate estimates 

B{B+ J/^K*+) = B{B+ ^ J/i;K+)Rf+, (25) 
B{B+ ^ J/i,K*+) = BiB"^ ^ J/^K') (^^^ <r, (26) 

BiB+ ^ J/V^ir*+) = BiB' J/i,K*') 1^-^^ (27) 
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We use for the first factor on the right-hand side of these estimates the world average values 
for the branching fractions and form the weighted average of these three measurements, 
thereby reducing the net statistical uncertainty. Because we employ in this calculation the 
world averages that have been determined assuming that = fd, these results depend 
implicitly on this assumption. Note that this estimate of B{B^ J/ipK*^) is statistically 
independent of the world average value for this branching fraction. 
Using this procedure, we obtain the branching fractions 





J/ipK+) 


= (0.82 ±0.18 ±0.07) X 10-^ 


(28) 


- 


J/ijK^) 


= (1.14 ±0.27 ±0.09) X 10"^ 


(29) 


B{B+ - 


> J/tljK*+) 


= (1.73 ±0.55 ±0.15) X 10-^ 


(30) 


BiB"^ - 




= (1.39 ±0.32 ±0.11) X 10^1 


(31) 



The statistical and systematic uncertainties have been estimated by weighting the relative 
contributions in the world average values and our data. The uncertainties in the world 
average branching fractions used in this calculation are dominated by the most recent mea- 
surements by the CLEO collaboration These uncertainties are limited by the size of 
the CLEO sample, and are therefore largely statistical and independent. We have exam- 
ined the stability of these estimates to different assumptions concerning the independence 
of the quoted systematic uncertainties. We find that our results and their estimated uncer- 
tainties are insensitive to possible correlations in the systematic uncertainties in the CLEO 
measurements. 



3. Comparison with Theory 

We have compared our measured ratios of branching fractions times fragmentation frac- 
tions with a calculation of the two-body nonleptonic decay rates of B mesons, using the 
factorization hypothesis, chiral and heavy quark symmetries and data from semileptonic D 
meson decays 0. We adjusted the predicted ratios by the world average B meson lifetimes 
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13| to correct for the observed lifetime differences of these three states. Although several 



recent theoretical calculations of these branching fractions exist, we have selected Ref. 
for this comparison as it predicts all the branching fractions for the five decays studied 
here. The other model calculations have been made with varying theoretical assumptions 
and observational constraints, but they generally predict ratios of branching fractions that 
are in reasonable agreement with each other and our observations. It should be noted that 
these calculations generally assume the validity of factorization as applied to nonleptonic B 
meson decays but they differ in many details, such as the magnitude and shape of the form 
factors for B meson decay and the experimental constraints employed in the calculations. 
In Ref. 0, the form factors are normalized to D meson semileptonic decay data and are 
assumed to be consistent with simple pole dominance. This assumption has been criticised 
recently [0,^ in the light of data on the observed polarization in the decay B J /ipK*. 

The results of the comparisons are shown in Fig. ^. In order to compare ratios involving 
S° decays, we have assumed /„ = fd and taken fg = (0.40 ±0.06) /„. The predictions agree 
well with the observed ratios of branching fractions for all the decay modes. 



B. Ratios of Meson Fragmentation Fractions 

The b quark fragmentation fractions have not been directly measured in a hadron collider 
environment. Our data allow us to constrain the ratios of these fractions. However, we note 
that the measured branching fractions of the B^ and B^ mesons have been determined 
assuming equal fragmentation fractions of b quarks produced in T(45') decays. We therefore 
cannot employ the world average values for these quantities in estimating /„, fa or fg. 
Instead, we will make specific assumptions concerning the branching fractions. 

Expressed in terms of fragmentation and branching fractions, the two ratios and 
R^l+ give the relations 

U - ""^^ BiBO J/^pKO) ' ^^^^ 
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U ~ - J/^ir*o) ■ ^^^^ 

Under the assumption that the ratios of branching fractions on the right-hand side of Eqs. ^ 
and ^ are unity (which is the result of most quark model predictions), the weighted average 



of these two quantities gives 

^ = 0.99 ±0.19 ±0.08. (34) 

This result is consistent with the hypothesis that b quarks hadronize equally often into 
and mesons. 

The strange meson fragmentation fraction fg is constrained by the ratios that involve 
the B^ final state, yielding the relationships 

A""^^- B{BO^J/^^) ' ^^^^ 

U ~ ^(5° - J/^0) ' ^^^^ 

A - ""^^ B{BO J/^0) ' ^^'^ 
/. _ , B{B^ ^ J/^K°) 

/. " "^^^ B{BO J/^0) • ^^^^ 
If we assume equal B~^ and 1?° decay rates to the J/ipK* final states, equal B^ and 1?° 
decay rates to the J/ipK final states and /„ = /d, we obtain the ratios 

= (0.24 ± 0.07 ± 0.02) ^'f,^2 (39) 

(/n, fd) B{B^g J/llj(f)) 

-(0.39 ±0.11 ±0.04) §1,^41^. (40) 



ifu, fd) ' ' ■ ■ ' i3(5° - J/V^0) ■ 
The probability of B^ meson production inferred from these data depends on the ratios 
of branching fractions in Eq. ( P5| ) and (0). We take for the ratios of these fractions the 
values predicted in ^ and correct for lifetime differences as discussed earlier. We find that 

f 

= 0.34 ±0.10 ±0.03. (41) 

[Jui Jd) 

In phenomeno logical fragmentation models, the probabilities /„, fd and fs are related to 
the relative probabilities of producing a uu, dd and ss quark pair in the quark fragmentation 
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process Measurements of the relative probabilities of strange meson to light meson 



production in e~^e~ and hadron-hadron collisions [jTT| and in deep inelastic scattering have 
yielded values in the range of 0.3 to 0.4. A recent compilation of these data has yielded a 
value of 0.29 ± 0.015 for the relative rate of strange quark production to up or down quark 
production |^|, which agrees well with our values measured in b quark fragmentation. Taken 



together, these measurements indicate that the rate of ss suppression in quark fragmentation 
is largely independent of energy and flavor of the quark initiating the fragmentation process. 

C. Fragmentation Fractions of B Hadrons 

The hadronization of b quarks produces both B mesons and baryons. Most B hadron 
decay models predict that virtually all B baryons produced during the fragmentation process 
will subsequently decay via modes that include a A° baryon, and so we make that assumption 
here. A measurement of the rate of A^ production in b quark fragmentation gives us a direct 
measure of f\^, the probability that a bottom quark will hadronize such that a A° baryon is 
produced. This, combined with our measurements of the ratios of fragmentation fractions, 
allows us to make a determination of the values of /„, fa and fs- 

Studies of A+ production in semileptonic b quark decays Il26|-p8| have yielded measure- 
ments of the product 

UMK^K^-'^iX). (42) 

The naive spectator quark model would predict that the inclusive A° semileptonic branching 



fraction to a final state with a charm hadron Xc is in the range of 0. 10. Reference suggests 
a possible range of 0.10-0.13. We have therefore chosen to use the inclusive branching fraction 
i3(A° Xci^ue) = 0.115, the central value of the theoretical prediction, to estimate f\^. 
Assuming this value and the three most recent measurements of A° production in b quark 
fragmentation, we find that = 0.096 ± 0.017. It is of interest to note that the observed 
ratio of meson to baryon production in minimum bias pp collisions at ^/s = 630 GeV, 
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6.4±1.1, yields a baryon production fraction of 0.14±0.02. This is consistent with the value 
determined from A° semileptonic decays even though these two fractions are not necessarily 
expected to be equal. 

Using the condition that the fragmentation fractions should sum to unity, assuming that 
the fraction of charm B hadrons fc<^i and all B baryons decay via A° intermediate states, 
then 

fu + fd + fs + fA, = l. (43) 

This can be rearranged to determine values for /„, fd and fs using /a^ and our measured 
ratios of fragmentation fractions: 

fu fu 
Jd Jd 

= ■ (46) 

Is Js 

We find that = 0.39 ± 0.04 ± 0.04, fa = 0.38 ± 0.04 ± 0.04 and fs = 0.13 ± 0.03 ± 0.01. 
These values are all proportional to the term (1 — /a J and are therefore relatively insensitive 
to our assumption concerning the semileptonic branching fraction. 



VI. CONCLUSION 

We have measured the ratios of branching fractions times fragmentation fractions for the 
five decay modes J/ijK+, B+ J/V^iT* (892) + , B^ J/^jK^, B^ J /i/jK* (892)° 

and 5° ^ 7/^^0(1020). 

We have used these measurements, with the assumption that /„ = fd and with fs — 
(0.40 ± 0.06)/„, to determine the relative branching fractions of the B mesons into the 
observed final states. We have made the first measurement of a B^ branching fraction to a 
final state with a J/ip meson, yielding 
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^ J/i^(t)) = (0.93 ± 0.28 ± 0.10 ± 0.14) x 10"^ 

We have also used our data in conjunction with the current world average branching fractions 
to find 



B{B+ - 




= (0.82 ±0.18 ±0.07) X 10"^ 


8(3" - 




= (1.14 ±0.27 ±0.09) X 10"^ 


B{B+ - 


. J/^IjK*+) 


= (1.73 ±0.55 ±0.15) X 10-^ 


B{B'' - 




= (1.39 ±0.32 ±0.11) X 10~^ 



These data are consistent with a B^ branching fraction approximately equal to those of the 
B~^ and B'^ decays into topologically similar final states. The observed branching fractions 
are in good agreement with model calculations employing factorization, chiral symmetry 
and heavy quark symmetries. 

An analysis of the ratios of branching fractions supports the widely-held assumption 
that the probabilities of producing 5+ and B^ mesons in b quark fragmentation are equal. 
We measured the ratio of fa to /„ to be 0.99 ± 0.19 ± 0.08. If we assume the theoretically 
predicted ratios of branching fractions for the B^ — > J/V^0(1O2O) relative to topologically 
similar decay modes for the S+ and S° mesons and we assume fu — fd: we determine 
that fs/fu — 0.34 ± 0.10 ± 0.03. Employing an estimate for the fraction of A° production 
in b quark fragmentation, f^^ = 0.096 ± 0.017, we determine /„ = 0.39 ± 0.04 ± 0.04, 
fd = 0.38 ± 0.04 ± 0.04 and f, = 0.13 ± 0.03 ± 0.01. Thus, our results imply a suppression 
of ss production relative to uu and dd production in b quark fragmentation similar to that 
measured in e"'"e~ and deep inelastic scattering experiments. 

VII. ACKNOWLEDGEMENTS 

We thank the Fermilab staff and the technical staff at the participating institutions for 
their essential contributions to this research. This work is supported by the U. S. Depart- 
ment of Energy and the National Science Foundation; the Natural Sciences and Engineering 

28 



Research Council of Canada; the Istituto Nazionale di Fisica Nucleare of Italy; the Ministry 
of Education, Science and Culture of Japan; the National Science Council of the Republic 
of China; and the A. P. Sloan Foundation. 



29 



REFERENCES 



[1] S. L. Glashow, Nucl. Phys. 22, 579 (1961); 
S. Weinberg, Phys. Rev. Lett. 19, 1264 (1967); 

A. Salam, in Elementary Particle Theory, edited by N. Svartholm (Almquist and Wik- 
sells, Stockholm, 1968), p. 367. 

[2] M. Wirbel, B. Stech and M. Bauer, Z. Phys. C 29, 637 (1985); 
M. Bauer, B. Stech and M. Wirbel, Z. Phys. C 34, 103 (1987). 

[3] A. Deandrea, N. Di Bartolomeo, R. Gatto and G. NarduUi, Phys. Lett. B 318, 549 
(1993). 

[4] G. Kramer and W. F. Palmer, Phys. Rev. D 46, 2969 (1992); 

G. Kramer and W. F. Palmer, Phys. Rev. D 46, 3197 (1992). 

[5] M. S. Alam et al. (CLEO Collaboration), Phys. Rev. D 50, 43 (1994). 

[6] H. Albrecht et al. (ARGUS Collaboration), Z. Phys. C. 48, 543 (1990). 

[7] F. Abe et al. (CDF Collaboration), Phys. Rev. D. 53, 1051 (1996). 

[8] S. Abachi et al. (DO Collaboration), Phys. Rev. Lett. 74, 3548 (1995). 

[9] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 75, 1451 (1995). 

[10] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 76, 2015 (1996). 

[11] V. Ammosov et a/., Phys. Lett. 93B, 210 (1980); 

R. Brandelik et al. (TASSO Collaboration), Phys. Lett. 117B, 135 (1982); 

P. K. Malhotra and R. Orava, Z. Phys.C 17, 85 (1983); 

W. Bartel et al. (JADE Collaboration), Z. Phys. C. 20, 187 (1983); 

H. Aihara et al, Phys. Rev. Lett. 53, 2378 (1984); 

G. J. Alner et al. (UA5 Collaboration), Nucl. Phys. B258, 505 (1985); 
R. E. Ansorge et al. (UA5 Collaboration), Z. Phys. C 41, 179 (1988); 



D. Bortoletto et al, Phys. Rev. D 37, 1719 (1988); 
F. Abe et ai, Phys. Rev. D 40, 3791 (1989). 

[12] P. Abreu et al. (DELPHI Collaboration), Z. Phys. C. 61, 407 (1994). 

[13] L. Montanet et al. (Particle Data Group), Phys. Rev. D 50, 1173 (1994). 

[14] F. Abe et al., Nucl. Instrum. Methods Phys. Res. A 271, 387 (1988); 
D. Amidei et al, Nucl. Instrum. Methods Phys. Res. A 350, 73 (1994); 
F. Abe et al. (CDF Collaboration), Phys. Rev. D 52, 4784 (1995). 

[15] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B327, 49 (1989); 

M. L. Mangano, P. Nason and G. Ridolfi, Nucl. Phys. B373, 295 (1992). 

[16] A. D. Martin, W. J. Stirling and R. G. Roberts, Phys. Rev. D 47, 867 (1993). 

[17] C. Peterson, D. Schlatter, I. Schmitt and P. M. Zerwas, Phys. Rev. D 27, 105 (1983); 
J. Chrin, Z. Phys. C 36, 163 (1987). 

[18] P. Avery, K. Read and G. Trahern, Cornell Internal Note CSN-212, 1985 (unpublished). 

[19] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 75, 3068 (1995). 

[20] T. Sjostrand, Comput. Phys. Commun. 82 , 74 (1994). PYTHIA Version 5.7 was used. 

[21] P. Abreu et al. (DELPHI Collaboration), Z. Phys. C 68, 353 (1995); 
M. Acciarri et al. (L3 Collaboration), Phys. Lett. B 345, 589 (1995); 
D. Buskulic et al. (ALEPH Collaboration), Z. Phys. C 69, 393 (1996). 

[22] M. Gourdin, A. N. Kamal and X. Y. Pham, Phys. Rev. Lett. 73, 3355 (1994). 

[23] R. Aleksan, A. Le Yaouanc, L. Oliver, O. Pene and J.-C. Raynal, Phys. Rev. D 51, 
6235 (1995). 

[24] G. Bocquet et al, Phys. Lett. B 366, 447 (1996). 

[25] V. V. Anisovich and V. M. Shekhter, Nucl. Phys. B55, 455 (1973); 



31 



R. D. Field and R. R Feynman, Nucl. Rhys. B136, 1 (1978); 

B. Andersson, G. Gustafson and B. Soderberg, Z. Rhys. C 20, 317 (1983). 

[26] R. Abreu et al. (DELRHI Collaboration), Z. Rhys. C 68, 375 (1995). 

[27] D. Buskulic et al. (ALERH Collaboration), Rhys. Lett. B 357, 685 (1995). 

[28] R. Akers et al. (ORAL Collaboration), Z. Rhys. C 69, 195 (1996). 

[29] T. Mannel and G. A. Schuler, Rhys. Lett. B 279, 194 (1992). 



32 



FIGURES 



10000 



CN 



7500 - 



> 
0) 



Ln 5000 - 



c 

> 



2500 - 




2.90 



3.00 3.10 

Gss lii'^/Li'' (GeV/c^) 



3.20 
2' 



3.30 



FIG. 1. The dimuon invariant mass distribution for the inclusive J/V' event sample. The shaded 
distribution is for same-sign dimuon candidates. 
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FIG. 2. The J/ijjK'^ invariant mass distribution is shown in (a). The invariant mass distribu- 
tions for the J/jpK^, J/tpK*^, J/ipK*^ and J/ipcp candidate events are shown in (b) through (e), 
respectively. 
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FIG. 3. The tt'^'tt" invariant mass distribution is shown in (a) for candidates. The K'^K~ 
invariant mass distribution is shown in (b) . The KgTr~^ invariant mass distribution is shown in (c) . 
The K'^Tr~ invariant mass distribution is shown in (d). 
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FIG. 4. The Kgir'^ invariant mass distributions for the S+ signal and sideband regions are 
shown in (a). The signal region is represented by the unshaded histogram. The sideband 
region, normalized to the non-5+ background in the 5+ signal region, is the shaded distribution. 
The corresponding distributions for the K'^tt" and K'^K~ invariant mass are shown in (b) and 
(c), respectively. 
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FIG. 5. Comparison of observed ratios of branching fractions (squares) with the theoretical 
prediction (triangles) described in the text. The error bars reflect the statistical and systematic 
uncertainties of the observed ratios added in quadrature. 
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TABLES 

TABLE I. The definition of the signal regions, sideband regions and the number of observed 
events associ atcd with the decays B+ J/ipK*+, J/ipK*^ and B^^^ J/i]j(l). The number 

of observed events is calculated in three ways: fitting the resonant structure in the Kit and KK 
invariant mass distributions either using the sideband distributions to estimate the background 
{Ngb) or using a second order polynomial to estimate the background (Nfu), and using the number 
of observed B candidate events, correcting for the loss of resonant decays due to the two-body mass 
requirement (A„,j„). The number of events obtained using the sideband subtracted background, 
Ngb, is used to calculate the branching fraction ratios. 

J/V'i^°7r+ J/i)K+TT- J/^PK+K- 

Signal Region (GeV/c^) 5.235-5.325 5.235-5.325 5.320-5.410 

Sideband Regions (GeV/c^) 5.000-5.220 5.000-5.180 5.100-5.305 

5.340-5.600 5.380-5.600 5.425-5.700 

(events) 21.3 ±6.1 119 ± 20 26.7 ±7.3 

Nfit (events) 17.0 ± 6.5 108 ± 27 27.3 ± 7.4 

Nyjin (events) 16.0 ± 5.3 119 ± 18 34.4 ± 7.3 
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TABLE II. The number of observed signal events and the reconstruction efficiencies for the five 
decay modes. The geometrical efficiency includes the meson Pt requirements and the acceptance of 
the tracking fiducial volume. It does not include the J/ip and light quark meson branching fractions. 
The efficiencies of the B proper decay length requirement and the reconstruction efficiencies of the 
light quark mesons are also listed. Some of the systematic uncertainties are correlated as they have 
common sources. These correlations are taken into account when ratios of the observed decay rates 
are determined. 



Channel 


Events 


^geom 






^meson 


B+ J/i)K+ 


154 ± 19 


(10.1 ±0.8) X 10" 


-2 


0.900 ±0.005 


0.979 ± 0.031 




36.9 ± 7.3 


(7.95 ± 0.58) X 10" 


-2 


0.876 ± 0.007 


0.857 ±0.012 


B+ J/'4;K*+ 


21.3 ±6.1 


(4.17 ±0.30) X 10- 


-2 


0.880 ± 0.006 


0.839 ± 0.015 




119 ±20 


(8.39 ± 0.59) X 10" 


-2 


0.893 ±0.006 


0.958 ± 0.032 




26.7 ±7.3 


(10.7 ±0.8) X 10" 


■2 


0.884 ± 0.020 


0.904 ± 0.058 



TABLE III. The systematic uncertainties in the relative efficiencies for the different channels. 

Effect Systematic Uncertainty (%) 
B meson confidence level requirements 1 
Detector simulation 5 
K'^ reconstruction efficiency 1 
Vector meson polarization 2.5 
reconstruction 1 
B Pt spectrum 1-5 
Effects of excited B meson production 1-4 
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TABLE IV. The ratios of fragmentation fractions times branching fractions for the various B 
meson final states. The ratio Rj is located in the ith row and jth column. The uncertainties are 

statistical and systematic, respectively. 

(■//V'i^°) (j/v>i^*°) jjm 

iJ/ipK+) 1.15 ±0.27 ±0.09 1.92 ±0.60 ±0.17 1.59 ± 0.33 ± 0.12 0.41 ± 0.12 ± 0.04 

{J/ipK^) 1.68 ±0.58 ±0.11 1.39 ±0.36 ±0.10 0.35 ± 0.12 ± 0.03 

{J/il;K*+) 0.83 ± 0.27 ± 0.07 0.21 ± 0.08 ± 0.02 

(J/ipK*^) 0.26 ± 0.08 ± 0.02 
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